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ABSTRACT 

We present the first weak-lensing-based scaling relation between galaxy cluster mass, Mwl, and 
integrated Compton parameter Kjph- Observations of 18 galaxy clusters at z :^ 0.2 were obtained 
with the Subaru 8.2-m telescope and the Sunyaev-Zel'dovich Array. The Mwl— Ksph scaling relations, 
measured at A=500, 1000, and 2500 pc, are consistent in slope and normalization with previous results 
derived under the assumption of hydrostatic equilibrium (HSE). We find an intrinsic scatter in AfwL at 
fixed Kjph of 20%, larger than both previous measurements of A/hse — i^sph scatter as well as the scatter 
in true mass at fixed l^ph found in simulations. Moreover, the scatter in our lensing-based scaling 
relations is morphology dependent, with 40% larger Mwl for undisturbed compared to disturbed 
clusters at the same Ygph at r^oo- Further examination suggests that the segregation primarily results 
from the inability of our spherical lens models to faithfully describe the three-dimensional structure 
of the clusters, in particular, the structure along the line-of-sight. We find that the ellipticity of the 
brightest cluster galaxy, a proxy for halo orientation, correlates well with the offset in mass from the 
mean scaling relation. This confirms empirically that line-of-sight projection effects are an important 
systematic uncertainty in lensing-based scaling relations. 
Subject headings: galaxies: clusters: general — gravitational lensing — cosmology: observations 
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Galaxy clusters are tracers of the highest peaks in the 
matter density field, and as such, their abundance as 
a function of mass and redshift depends strongly upon 
cosmology. Along with data from a variety of other tech 



niques (Komatsu ct al. 2011; Hicken et al.||2009 Kessler 



et al.|[2009[ jRercival et al.||2010| ) , galaxy cJ 

have placed precise constr aints on the ACDM cosmo- 



uster surveys 



logical parameters f Vikhlininetalj[2009b ; Mantzetal. 
2010b Vanderlinde et al.,201(j|. In order to improve the 



constraints from cluster surveys, the scaling relationship 
between survey observables and the total cluster mass 
must be precisely determined. 

Surveys at millimeter wavelengths can be used to iden- 
tify large numbers of clusters via the Sunyaev-Zel'dovich 
(SZ) effect, as ha s been demonstrated by the South Pole 
Telescope (SPT; 'Staniszcwski et al.||2009| IVanderlinde 



et al.||2010| [Wi lfiamson ct al. 20 ^7 the Atac ama Cos 



mology I'elescopc (ACL; Marriage et al. '2011), and the 
Planck satellite (Planck Collaboration et al.^^2011} . The 
SZ effect is a spectral distortion ot the cosmic microwave 
background (CMB) resulting from the Compton scatter- 
ing of CMB photons by the hot electrons of the intra- 
cluster medium (ICM). The magnitude of this distortion 
is determined by the integral of the pressure along the 
line of sight through the cluster. The volume integral of 
pressure is the total thermal energy content of the ICM, 
whic h should directly trace the depth of the potential 



well (Carlstrom et al. 2002), so a tight scaling between 



the SZ signal (Y , see Section 3.2 ) and mass is expected 



A low-scatter correlation is indeed found in simulations 
of galaxy clusters 



(e-g-, 



da Silva et al. 2004 Motl et al 



2005; Nagai 2006| |StaneJt et al-pOlOf ; aHd is suggested 



in the best comparisons ot SZ o bservations and cluster 
mass estimates available to date (Bonamente et al.|2008 



Marrone et al. 



Plagge et al.[2010 



Andersson et al.|2011 ). However, the 



mass-y scaling relation remains poorly determined from 
these observations, and the uncertainty in the relation 
is already limiting the cosniolo gical constraints derive d 
from small SZ-selected samples (Vanderlinde et al. 2010). 
Moreover, previous comparisons ot Y with mass have 
relied on mass estimates derived assuming hydrostatic 
equilibrium (HSE). This assumption is certain to result 
in at least some bias, although the degree of bias is poorly 
known; observational constraints suggest ~ — 20% 
within radii of interest (|Mahdavi et alj2008 



Zhang et al. 



arger biases 



2010 1 , though some simula tions have found 
(e.g., Cavalier e et al.|[20TT |. 

One promising technique for improving our knowledge 
of the mass-y rela t ion is weak gravitational lensing (e.g 



Dahle et al. 2002 Hoekstra 2007 Okabe et al. 2010a) 



Unlike other techniques, weak lensing directly probes the 
gravitational potential of the cluster, providing a way 
to test for mass biases in other mass estimators. How- 
ever, the weak lensing signal is sensitive to the total mass 
projected along the line-of-sight, and thus likely suffers 
projection-related errors in mass measurements. Several 
authors have calibrated X-ray observables against lens- 
ing mass estimates, typically finding consistency with 
(but with larger scatter than) X-ray -only scaling rela- 



tions (Smith et al. 



2005 



Bardeau et al. 2007; Hoekstra 
20071 | /^liang et al."! S0071l200^1|Rykotf et al.||2008[ |Ukabel 
et al.| [2010b Hoekstra et al.||2011|. This technique has 



not yet been used to calibrate the SZ observable except 
in the cluster core (Marrone et al.^^2009J . An important 
next step is to make such a comparison at the larger 
radii typically used in cosmological studies with galaxy 
clusters. 

The calibration of mass-observable relations is one of 
the key goals of the Local Cluster Substructure Survey 
(LoCuSa^j), a multi- wavelength survey of galaxy clus- 
ters at 0T5 < z < 0.3 se lected fr om the J^Qg^T All-sk y 
Survey jEbehng et al.|19 98, 2000] [Bohringer et al.|2004| . 
The LoCuSS sample is ROSAT X-ray luminosity limited 
(Lx[0.1-2.4 keV] > 2 x 10"^"' ergs~^); the large dynamic 
range in Lx and lack of morphological selection ensure a 
morphologically diverse population for calibration stud- 
ies. Weak lensing mass estimates are available for an 
initial sar nple of 29 LoCuSS clusters ( Okabe et al.|2010a 
hereafter OklO). In this paper, we combine these esti 



mates with SZ data from the Sunyaev-Zel'dovich Array 
(SZA) in order to constrain the mass-F scaling relation 
and its scatter. We also explore the possibility of bi- 
ases in the mass estimates related to cluster morphology, 
which our data suggest is an important systematic. 

Section [2] describes the sample of clusters as well as 
the weak lensing and SZ observations. In Section [Sj we 
discuss the analysis techniques used to derive the weak 
lensing mass estimates and determine the SZ observables. 
Our scaling relation results are presented in Section |4] 
and our interpretation of the findings is discussed in Sec- 
tion [5l Our conclusions are reviewed in Section |6l We 
assume a flat ACDM cosmology: J^m = 0.27, Ha = 0.73, 
Ho = 73 km s~^ Mpc^^. 

2. OBSERVATIONS 



2.1. Cluster Sample 

The clusters considered in this work are drawn from the 
LoCuSS sample (G. P. Smith, in prep.), which was se- 
lected from the ROSAT All-Sky Survey according to the 
redshift and Lx cuts above, within the declination range 
-70° <6< -h70° and with uh < 7 x lO^^cm-^. A first 
batch of "LoCuSS clusters" satisfying this selection func- 
tion, and visible from Mauna Kea on the assigned nights, 
was observed with Subaru/Suprime-Cam. Of this batch, 
21 clusters meet the additional criteria that the V- and 
i'-band data is of the quality required for reliable weak 
lensing mass measur ements, and that a single NFW halo 
ields an acceptable fit to the data 



(Navarro et al. 



19961) 



( OklO[ see Section |3.1[ ). These 21 clusters make up our 
sample (Table IlJ . I'he sample spans a factor of 4 in X- 
ray luminosity and includes clusters with a broad range of 
X-ray morphologies, includi ng cool-core (hereafter CC) 



clusters (e.g., A 383, A 1835: Smith et al. 2001 Schmidt 



et al. 2001 ), cold fron t clusters (e.g., RXJ 1720 1-^2638: 
Mazzotta et al. 2001), and clusters known to be under- 
going mergers (e.g., A 209: Giovannini et al. |2009j and 
references therein). 

2.2. Gravitational Lensing 

The acquisition and processing of the Subaru/Suprime- 
Cairp^data upon which the weak lensing mass measure- 
ments are based are described in detail by OklO Briefly, 
the data were acquired in the V- and i'-bands in excellent 
atmospheric conditions FWHM ~ 0.7", consistent with 
the stringent requirements for precise measurement of 
galaxy shapes. Weakly-lensed background galaxies were 
selected to have a minimum color offset from the clus- 
ter red sequence in the color- mag nitude plane fo llowing 
Umetsu fc Broadhurst] (|2008|) and [Umetsu et al] (|2009 ) 



Background galaxy shapes were measured using KSB 
( jKaiser et al.,1995) , achieving a residual mean ellipticity 
of point sources after removal of the point spread func- 
tion of ~ 10""', and a typical number density of galaxies 
for use in the weak lensing analysis of ~ 10—30 arcmin"^. 

2.3. Sunyaev-ZeVdovich Effect 

The sample of 21 clusters was observed with the SZA, 
an 8-element radio interferometer optimized for measure- 
ments of the SZ effect. The SZA is a subset of the 
23-element Combined Array for Research in Millimeter- 
wave Astronomy (CARMA). Data were obtained in an 
8 GHz passband centered at 31 GHz between 2006 and 
2009 (Table fTl) . Over the period of these observations, 
the array was sited at two different locations and oc cu- 
pied three different c onfigurations, as described in |CuF] 



verhouse et al. 



http : //www . sr . bham . ac . uk/locuss 



( |2010| ). For most of this period, the SZX 

was configured with six elements in a compact array 
sensitive to arcminute-scale SZ signals, and two outrig- 
ger elements providing discrimination for compact radio 
source emission. The compact array and outrigger base- 
lines provide baseline lengths, in units of number of wave- 
lengths A, of 0.35—1.3 kA and 2—7.5 kA, respectively. 
Two of the clusters (Abell 209 and ZwCl 1454.8-1-2233) 
were also briefly observed with an eight-compact-element 

^^ Based in part on data collected at Subaru Telescope and ob- 
tained from the SMOKA, which is operated by the Astronomy 
Data Center, National Astronomical Observatory of Japan. 
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array configuration, with greater SZ sensitivity but little 
intrinsic power to distinguish SZ signal from contamina- 
tion. Integration times were tailored to the magnitude 
of the SZ signal in each cluster and, secondarily, to the 
degree of contamination from radio sources. The total 
on-source unflagged integration times range from 5 to 
69 hours, and the corresponding noise levels achieved by 
the short (SZ-sensitive) baselines were 0.14 to 0.45 mJy 
(12 to 40 ^K), with synthesized beams of ^2' FWHM. 
Absolute calibration of the data is determined through 
periodic me asurem ents of Mars, calibrated against the 
Rudy et al. (1987) model for this source. A systematic 
uncertamty ot 57o applies to the SZ measurements due 
to the uncertainty in the absolute calibration, this is not 
included in the reported errors. 

Exclusion of emissive sources from the SZA data is 
most simply achieved when the sources are spatially com- 
pact. For such sources the emission measured on the 
longer outrigger baselines (smaller angular scales) can be 
assumed to be representative of the emission measured on 
the short, SZ-sensitive baselines (larger angular scales) 
and thereby separated from the SZ signal. Sources that 
are not point-like on the scale of 20" leave residual emis- 
sion in the short baseline maps that fills in the cluster SZ 
signal. Three clusters in the sample suffer from extended 
radio source contamination: 

• Abell 115 hosts the bright radio source 3C 28, 
which is >100 mJy at 5 GHz and has lobes ex- 
tending more than an a rcminute across the cluster 
(Giovannini et al. 1987). While there is some dis- 
crimination on the short baselines between the ex- 
tended central source and the more extended SZ 
decrement. Reliable separation of these two re- 
solved structures is not possible in these data. 

• RX J1720. 1+2638 contains a bright source near the 
cluster center that is resolved on sc ales of 30" at 
1.4 G Hz in the VLA FIRST survey ([Becke r et al 
1995). The SZA data suggest that this source is 



also extended at 31 GHz, where the SZ decrement 
appears to be split through its center by an emis- 
sion source after a point-like source model is deter- 
mined from the long baselines and removed. 

• Abell 291 hosts three 31 GHz sources within an ar- 
cminute of the pointing center. One of these is re- 
vealed to be a double-lobed radio source in FIRST, 
and residual emission at this position appears to fill 
in the SZ decrement, leading to an offset of >20" 
between the SZ centroid and BCG position. 

Observations of these objects in the 90 GHz SZA band 
would likely provide greater spectral discrimination be- 
tween the radio source emission and SZ signal, this possi- 
bility will be examined through future observations. For 
this work, we omit these three clusters from further anal- 
ysis, leaving 18 clusters with which to study the mass-K 
scaling relation. 

3. ANALYSIS 
3.1. Weak Gravitational Lensing 
The procedures for determining galaxy cluster masses 



and related radii are described in ful l by|OklO 



galaxy sample - i.e., we use only those galaxies that are 
redder or bluer than the cluster red sequence. Exclusion 
of unlensed galaxies in this manner is essential because 
cluster/foreground galaxies can bias A/2500 and Afsoo low 
by ^ 50 — 20%. The mean redshift of the color-selected 
background galaxies is estimated by matching their col- 
ors and magn itudes to the COS MOS photometric red- 
shift catalog (Ilbert et al. 20091. The tangential shear 
signal is then measured in logarithmically spaced bins, 
the measurement in each bin being the error-weighted 
mean tangential shear of the galaxies in that bin. The 
bins are centered on the brightest cluster galaxy (BCG) 
of each cluster. 

The spherical mass (or spherical overdensity mass) A/a 
is defined as the mass enclosed within a radius ta for 
which the average interior density is A times the criti- 
cal density of the universe {pc = 3//(z)^/87rG). A/a is 
estimated by fitting the measured shear profile to the 
NEW model prediction parameterized by A/a and con- 
centration CA, where the NEW mass profile ( [Navarro 
1996) is given as p ex r^^(l + ca^'/^a)"^ Both 
and CA are allowed to be free in the fit — no mass- 
concentration rela tion is assumed , as is often done in 
other works (e.g., Hoekstra 2007) — and the A/a values 
for each cluster are determined by marginalizing over 
the posterior distribution of concentration values. The 
spherical weak lensing masses at the three overdensities 
used in thi s work , A =2500, 1000, and 500, are listed in 
Table 8 of [OklO Overdensity radii, which are calcula- 
ble from the masses (A/a = Airr^Apc/S), are given in 



et al. 



Table [2l 



3.2. Integrated Sunyaev-Zel'dovich Effect Signal 

For scaling relation analyses, the SZ signal is typically 
quantified using the Compton j/-parameter integrated 
within a region centered on a cluster: 



Y= ydn^ 
1 ut 



1 kB<JT 



D 



rar.c^ 



dl 



oT.dA 



P{r)dA. 



(1) 



^ we sum- 
marize salient details here. We adopt [OklOf s "red-|-blue" 



The quantity YD\ is called the intrinsic y-parameter, as 
it removes the distance dependence of the right side of 
Equation [l] 

It is customary to integrate Compton-y over the solid 
angle subtended by the cluster, as the average change 
in the CMB temperature within an aperture (AT ex 
YTqmb) is expected to be a robust observable. One stan- 
dard method of accomplishing this is to compare the ob- 
served sky signal with a parameterized radial SZ profile, 
projected to two dimensions and filtered in the same way 
as the sky data. However, both in interferometric data 
such as ours and in data from single-aperture telescopes, 
the large-scale SZ signal is attenuated by spatial filter- 
ing and/or confuse d with the a nisotropy of the primary 
CMB (e.g., Plagge et al.|[20l"o| . The poor constraint on 
the large-scale behavior o± the profile is coupled into the 
derived Y as additional uncertainty from the unknown 
large I (line of sight distance) behavior of the profile. 

This effect can be partially mitigated by adopting as 
an observable the integral of the model profile over a 



Marrone et al. 



TABLE 1 
Cluster Sample and SZ Observational Parameters 



Cluster 


R.A.^ 


Dec'' 


z 


W 


tint 


e c 


Resolution'^''^ 


rp C 

J-rms 










(lO^-^ergs-i) 


(hrs) 


(mjy) 


(") 


{^^K) 


Abell 68 


00:37:05.90 


+09:09:26.6 


0.255 


8.8 


17.3 


0.22 


122x139 


17 


AbellllSN 


00:55:50.38 


+26:24:36.0 


0.197 


8.6 


7.7 


0.36 


126x130 


29 


Abell 209 


01:31:53.47 


-13:36:46.1 


0.206 


7.3 


19.3 


0.21 


91x125 


24 


RXJ0142.0+2131 


01:42:02.64 


+21:31:19.2 


0.280 


5.9 


11.6 


0.30 


114x126 


27 


Abell 267 


01:52:41.93 


+01:00:24.1 


0.230 


8.1 


18.4 


0.30 


115x160 


21 


Abell 291 


02:01:43.11 


-02:11:48.1 


0.196 


5.7 


27.8 


0.22 


101x124 


23 


Abell 383 


02:48:03.50 


-03:31:45.0 


0.188 


5.3 


25.0 


0.27 


128x159 


18 


Abell 521 


04:54:06.90 


-10:13:24.6 


0.247 


9.5 


12.8 


0.32 


117x210 


26 


Abell 586 


07:32:20.31 


+31:38:02.0 


0.171 


6.6 


13.1 


0.25 


120x137 


21 


Abell 611 


08:00:56.74 


+36:03:21.6 


0.288 


8.1 


30.2 


0.25 


120x130 


21 


Abell 697 


08:42:57.80 


+36:21:54.5 


0.282 


9.6 


15.8 


0.35 


118x128 


30 


Abell 1835 


14:01:02.02 


+02:52:41.7 


0.253 


22.8 


13.5 


0.36 


117x149 


27 


ZwCl 1454.8+2233 


14:57:15.09 


+22:20:34.2 


0.258 


7.8 


68.6 


0.14 


119x123 


12 


ZwCl 1459.4+4240 


15:01:23.13 


+42:20:39.6 


0.290 


6.7 


21.7 


0.22 


115x131 


19 


Abell 2219 


16:40:22.60 


+46:42:22.0 


0.228 


12.1 


27.3 


0.20 


114x133 


17 


RXJ1720.1+2638 


17:20:09.90 


+26:37:27.8 


0.164 


9.5 


26.9 


0.19 


108x125 


19 


Abell 2261 


17:22:27.08 


+32:07:58.6 


0.224 


10.8 


9.9 


0.31 


114x134 


26 


RXJ2129.6+0005 


21:29:39.90 


+00:05:19.8 


0.235 


11.0 


24.5 


0.27 


113x167 


19 


Abell 2390 


21:53:36.70 


+17:41:31.0 


0.233 


12.7 


4.8 


0.45 


109x136 


40 


Abell 2485 


22:48:31.13 


-16:06:25.6 


0.247 


5.9 


21.2 


0.19 


106x126 


18 


Abell 2631 


23:37:38.80 


+00:16:06.5 


0.278 


7.9 


16.1 


0.28 


140x152 


17 



Pointing center for SZA observations. 
Unflagged, on-source time. 
Short, SZ-sensitive baselines only. 



Synthesized beam size. 

spherical volume, 



-* sph — 



1 aT 



P{r)dV. 



(2) 



Fsph is relatively insensitive to the unconstrained modes 
in the data, and can be calculated vifithout resorting to 
arbitrary outer r adius cutoffs for the line-of-sight integra- 
tion (as noted by Arnaud et al.|2010 ) . Furthermore, such 
integration is standard in the analysis of X-ray measure- 
ments of clusters (e.g., Zhang et al.|2008" Vikhlinin et al. 
2009a). We therefore perform our scaling relation analy- 
sis using the spherically integrated intrinsic y-parameter, 

Since interferometers do not measure the total sky sig- 
nal, but instead sample the Fourier transform of the sky 
signal only at the spatial frequencies determined by their 
uv coverage, we must restore missing spatial information 
in order to determine our observable. We accomplish this 
by adopting a model profile to fit to our visibility data; 
specifically, the generalized NFW (G NFW) pressure pro- 
file proposed bylNagai et al. (2007), 



P{r) = 



Pn 



a;7(i + a:-)(^-^)/"^ 



(3) 



where x = r/rg. The model contains five parameters: 
three exponents (a, /3, 7) which determine the shape 
of the profile, a normalization Pq , and a scale radius Vg 
(alternatively, a concentration parameter). 

Our data are not capable of placing useful constraints 
on the three shape parameters due to degeneracies be- 
tween them, the moderate signal to noise, and the lim- 
ited range of angular scales probed by the interferometer. 
We fix the v alues of the shape parameters to the |Arnaud| 
et al I ( |2010[ ) average values, [a, ^, 7] = [1.05, 5.49, 0.31J, 
which were determined from simulations and scaled X- 
ray pressure profiles of galaxy clusters in the low-redshift 



(z < 0.2) REXCESS sample. We estimate the remaining 
pressure model parameters (Pq, Ts, and the centroid posi- 
tion) and their uncertai nties using a Markov cha in Monte 
Carlo method (MCMC; |Bonamente et al.|2006[ ). At each 
step of the MCMC, we integrate the spherical profile 
along the line of sight, filter it with the primary beam 
pattern of the SZA antennas (approximately a gaussian 
with FWHM 11.0'), and Fourier transform for compari- 
son with the measured Fourier components (visibilities). 
The positions and fiuxes of detected radio sources are in- 
cluded in the MCMC as additional free parameters, and 
are marginalized over to determine the pressure model 
parameters. 

The centroid of the SZ signal is a free parameter in 
our model to allow for the possibility of real offsets be- 
tween the center of the ICM and the BCG. The typ- 
ical uncertainty in the best-fit SZ centroid is O.OSrsoo, 
corresponding to ^10". For all 18 clusters, the best- 
fit SZ centroid is within 0.15r5oo of the BCG position 
(the assumed center of shear), with the SZ centroid in 
75% of our sample lying within 0.04r5oo of the respective 
BCGs. Larger offsets were typically found in clusters 
that lack a cool core, or show signs of merger activity. 
For example, the most significant offset is seen in A 2390 
('"Sz-WL = 0.15±0.04r5oo), which shows both a cool core 
and merger activity. After accounting for the measure- 
ment uncertainty, the distribution of offsets is reasonably 
consistent with recent X-ray /opt ical studies ( Sanderson 
(et al.,2009(|Haarsma et al.||2010[ ). 

lb calculate Ygph tor each cluster, we integrate the re- 
spective ICM model (Equations [2] &l3| over the spherical 
volume defined by the overdensit y ra dii obtained from 
the weak lensing models (Section |3.1[ Table l2| . These 
radii are not constrained from the SZ data alone, so 
a completely independent determination of Ygph at the 
overdensities of interest is not possible. The correlation 
between AfwL arid l^ph introduced by this procedure is 
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discussed in SectionjO] Note that the shared value of ta 
is the only coupling between the measurements of l^sph 
and the lensing analysis — there is no joint fitting of the 
mass and pressure profiles. The final Ygph measurements 
are listed in Table |2l 

3.3. Morphological Classification 

One of our aims is to check whether the shape, nor- 
malization, and intrinsic scatter of the mass-Fsph relation 
depends on cluster properties. Of particular interest is 
whether clusters might exhibit different behavior based 
on their dynamical state, which we often classify simply 
as merging or non-merging; however, dynamical state 
is difficult to estimate from SZ and weak lensing mea- 
surements alone. Here we employ a simple and broadly- 
accepted estimator of dynamical state: the centroid shift 
parameter (w) , defined as the standard deviation of the 
projected separation between the X-ray peak and the 
centroid of emission calculated in circular apertures cen- 
tered on the X-ray peak. 

We measure (w) for each cluster in our sample using 
archival Ch andra data, following the method of |Mohr| 



et al. ( 1995 1 as implemented by Maughan et al. 
i'he analysis is performed on exposure-corrected images 
with sources masked. The core (< 30 kpc) is excised 
from the calculation of the centroid, and the radii span 
0.05r5oo to rsoo in steps of 5% of r^oo- The core is in- 
cluded in the calculation of the X-ray peak. The error on 
(w) is derived from the standard deviation observed in 
cluster images resimula ted with Poisson noise, following 



Bohringer et al. 



(20101. Previous observational studies 



ot the (w) distribution (e.g., Maughan et al. 



2008 Pratt 



et al.|2009[|Bohringer et al.|20iup have otten divided clus- 
ters into "disturbed" and undisturbed" sub-samples at 
(w) = 10~^r5oo. We adopt this value to aid comparison 
of our results with the literature. Six of the 18 clusters 
are classified as disturbed {{w) > 10~^r5oo), and the re- 
maining 12 as undisturbed {{w) < lO^'^r^oo', Table [2||^ 
Note that although much of our discussion makes use o: 
this binary morphological classification system, equat- 
ing the disturbed systems with mergers and undisturbed 
systems with non-mergers, the true relationship between 
centroid shift and dynamical state is likely to be more 
complex. 

4. RESULTS 

The weak-lensing mass and Ygph measurements are 
compared in Figure [l] and a positive correlation is ev- 
ident at all radii. In this section, we define the scaling 
relation model that we will fit to the data, describe the 
regression analysis, and present the main results. 

4.1. Self-similar Scaling 

Without significant influence from non-gravitational 
processes, the mass, temperature, size, and other proper- 
ties of galaxy clus ters are expe cted to follow self-similar 
scaling relations ( [Kaiser 1986). These power-law rela- 
tionships follow directly troiri the assumptions of hy- 
drostatic equilibrium and an isothermal distribution of 

^^ Two of the three clusters excluded in Section |2.3| have X-ray 
data as well. Of these, one would qualify as disturbed (Abell 115N, 
{«)>=7.1xl0-2) and the other as undisturbed (RX J1720.1-I-2638, 
(ui}=0.25xl0~^), so their exclusion does not preferentially affect 
one subsample 



baryons and dark matter (e.g., Bryan fc Norman||1998 ), 
and form a useful reference point tor measurements of 
scaling relations. 

The simplest self-similar model, as derived in the 
above-cited papers, predicts that cluster mass {M) and 
temperature (T) are related as 



T ex M2/3£;(^)2/3aV3 



(4) 



where E{z) describes the redshift evolution of the Hub- 
ble parameter {E{z)^ = H{zf/Hl = f]M(l + zf + Vlp, 
for a flat cosmology). To derive the expected scahng 
between mass and the SZ signal, we note that the dou- 
ble integral on the right hand side of Equation [l] can 
be written as AfgasTo for an isothermal ICM, or cquiva- 
lently as fgasMT^ for /, 



gas 



M, 



with EquationHl we find the Mi 

M (X f-^J'A-'/^ 



gas/M. Combining this 
WL— ^sph scaling relation: 



YDiEiz 



-2/3 



3/5 



(5) 



This equation defines the self-similar reference model for 
our scaling analysis. 

4.2. Regression 

We analyze the scaling between mass and Ygph at three 
overdensity radii (A=500, 1000, and 2500) determined 
from the weak lensing measurements. We fit the data 
with a power law of the form 



M(rA) 
1014M, 



10-^ 







Y,pi,DlE{zry^ 

10-5Mpc^ 



(6) 



The regression is performed in linearized coordinates by 
using the base-10 logarithm of the data points. 

Our data have significant variations in their mea- 
surement precision, and simulations suggest that the 
M^i/h—Ysph correlation should have intrinsic scatter, 
placing significant demands on our data fitting tech- 
nique. The problem of linear regression with uncertainty 
in both axes and intrinsic scatter is complex, and several 
methods of parameter estimatio n have been formulated 
for this circumstance (e.g., Akrit as fc Bershady [1996 



Tremaine et al. 2002 Weiner et al. 2006). MethodsHiat 



Ignore the intrinsic scatter, particularly when accompa- 
nied by heterogeneous mea surement er rors, lead to bi- 
ased regression parameters ( Kelly|[2007 ) . 

We perform our linear regression using a three- 
parameter model: normalization A, slope B, and intrin- 
sic scatter aMiy- Parameter determinations are achieved 
through the Bayesian reg ression method presented in 
Kelly (2007). In that work Kelly demonstrates that this 
method performs well in regimes where parameter esti- 
mates from other methods (OLS, BCES, FITEXY) are 
significantly biased. The regression i s perfo rmed with 
publicly available IDL code written by [Kellyl 

4.3. Scaling Relations 

We fit the scaling relation (EquationJGl) to the full sam- 
ple of clusters with the slope _B as a free parameter, refer- 
ring to the result as the "free slope" (FS) fit (Table Isl. 
We also report results with the slope fixed to the self- 
similar value of i? = 0.6, which which we refer to as the 
"self-similar" (SS) fit. The posterior distributions of the 
normalization and scatter in the SS fit are obtained by 
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TABLE 2 
Cluster Mass and Kjpii 



Cluster (w) '■2500 -^^2500 1'sph,2500-DA '"1000 MiooO i'sph,100o£'A ''500 M500 i'sph.SOoi'A 

(10^^ rsoo) (Mpc) (10" Mq) (10^^ Mpc^) (Mpc) (lo" Mg) (10^^ Mpc^) (Mpc) (lo" Mq) (10"'^ Mpc^ 



Abcll 68 


1.17+0.14 


0.41 


, 07+0.57 


2.18tl\l 


0.70 


2.6912:™ 


^■5itl:ll 


Abcll 209 


0.68+0.10 


0.48 


2 11+°"^ 
^■^J^-0.45 


2.93l°:f3 


0.88 


5.05t°o:S 


em+_\:tl 


RXC J0142. 0+2131 


0.76+0.13 


0.45 


-'-■'''-0.30 


„ oc+0.27 
^•''°-0.27 


0.72 


2.95t°o:S 


4„o + 1.07 


Abell 267 


1.68+0.05 


0.42 


, „o + 0.24 

1-38-0.24 


1 67+°'^'^ 

J^-°'-0.15 


0.67 


2 31+°*1 


r, QO + 0.56 

2-88_o,5(, 


Abell 383 


0.20+0.04 


0.45 


1 fiq+0-23 
l.Dy_Q^23 


94+0.13 

'J-''*-0.12 


0.70 


„ ro+0.41 
Z.OO_Q 3g 


, ns+°-22 

l-'J8_o.l9 


Abell 521 


5.16+0.12 


0.38 


i.06iS:f, 


, r,„+0.12 


0.67 


^ o'7+0.46 
2.3^_0,46 


(,,+0.43 


Abell 586 


0.26+0.06 


0.57 


3.30tS:f, 


2 64+°™ 

^■°*-0.57 


0.89 


5-02ll:?^ 


3.72;;-^= 


Abell 611 


0.33+0.08 


0.44 


i-78tS:r, 


1 91+015 

^•"^-0.15 


0.75 


oQ+0.60 
■^■■^"-0.58 


3.66;°f, 


Abell 697 


0.58+0.09 


0.48 


2.i9tS:=? 


4.88lS:r, 


0.85 


4 96+°-^* 
*-»"-0.73 


13-43l^:fo 


Abell 1835 


0.22+0.01 


0.52 


„ »~+0.55 
^■' '-0.56 


5.92lS:|° 


0.91 


5 02+0-92 
''■-'^-0.87 


10.76+1-55 


ZwCl 1454.8+2233 


0.37+0.02 


0.35 


0-86_o:4o 


_„+o.os 
'J-' '-0.08 


0.60 


l-6»_o:52 


1 fii+0-2-* 

1^-01-0.23 


ZwCl 1459.4+4240 


1.29+0.16 


0.44 


1 Vo + 0.40 
-'-■'J_0.41 


2 48+°'2'' 


0.70 


on + 0-64 


.; 42 + 0.89 


Abcll 2219 


1.75+0.03 


0.58 


po + 0.57 

0.00_Q gQ 


D.^5_(, 45 


0.91 


5.83_(,;gj 


12.20i;-f2 


Abcll 2261 


0.81+0.07 


0.56 


42+0-42 
3-42_o 43 


„ o,+0.60 
2.81_(,42 


0.91 


5 72+0-71 

''-'^-0.68 


04+1.23 


RXC J2129. 6+0005 


0.43+0.11 


0.41 


1 00+0.51 
J^-33_o 52 


1 OO+0.13 


0.72 


2.86+°-^^ 


4 19+0 ■^^ 


Abell 2390 


0.78+0.03 


0.54 


3.03+S;S 


4 64+°'"' 


0.87 


5.0llg:™ 


8 29+^10 


Abcll 2485 


0.53+0.14 


0.37 


0.98+S:« 


i-i7l°o:l°o 


0.64 


2.04«:r, 


2 48+0-37 
^-*°-0.34 


Abcll 2631 


1.68+0.14 


0.49 


„ oQ + O-34 
^■■^•^-0.36 


2.36t°o:f3 


0.76 


3-58t°o:S 


n 2,+0.91 



1.01 
1.32 
0.99 
0.94 
0.96 
0.99 
1.22 
1.07 
1.26 
1.33 
0.86 
0.97 
1.27 
1.27 
1.05 
1.21 
0.93 
1.05 



,+1.18 
^-1.03 
3 + 1.26 

-1.15 
,+0.80 

-0.71 
" + 0.66 

-0.59 
1+0.66 

-0.57 
,+0.70 

-0.65 

n+1.89 

-1.54 
s + 0.91 
^-0.84 
4+1.18 
*-1.09 
a+1.64 

-1.46 
,+0.93 

-0.78 
^ + 0.97 

-0.86 
7+1.44 

-1.31 
-1 + 1.18 

-1.08 
-1 + 1.05 

-0.94 

6.811;-- 
-+0.86 

" -0.76 
- + 0.70 
^-0.67 



4.01] 
8.49I 

3.91J 
3.16J 
3.24J 
3.81^ 
6.50J 
4.98J 
8.041 
9.281 
2.51J 
3.761 
7.77I 
7.80I 
4. 50-* 



3.15: 
4.66^ 



6.57" 



+ 1.86 
-1.39 

io.i8ii:« 

c 00 + 2. 18 
O-JJ-i.ss 

„ + 1.10 



3.78^ 
1.151 
4.60I 
4.37-' 



0.8 



,; + 0.29 

-0.22 

n+1-23 

-0.94 

7+2.85 

-1.52 

5.06l}-^2 
23 lfi+7 00 

14 05+^-^" 

2-34lS:f4 
8-50l-iJ 

17.58ll:°I 
+ 1.S0 
1.13 
+2.01 
1.37 
+5.75 
3.55 

+ 0.87 

0.71 

+ 2.30 

-1.80 



4.40 
6.39 
11.15 
3.61 

8.37 



Note. — Cluster masses as reported in|QklO|fQr NFW halo fits to the shear profile. 
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Fig. 1. — The mass-y scaling relation at three overdensities, shown in the upper left corner of each plot. Two fits are shown at each 
overdensity: the best-fit line with a free slope is shown in black, and the best-fit line with the slope fixed to self-similar (l/B=5/3) is shown 
in red. The Ic uncertainty regions are shaded and marked with dashed lines . Results from prior analyses are plotted for co mparison in 
blue: the M— Kjpij scaling relation derived from the data oflBonamente et aLlpOOSJ (left panel), and the scaling relation from lAnderssonl 



|et al.| | |2011| (right panel) 

taking links with the appropriate slope from the markov 
chains generated by the fitting routine. 

At all three overdensity radii, the FS mass-l^ph rela- 
tions are slightly shallower than the self-similar predic- 
tion oi B = 0.6 (Figure [1]), although the differences from 
self-similarity are not statistically significant; the relation 
at A = 500 is the most discrepant at 1.40- significance. A 
jackknife test of the rsoo data, dropping one cluster at a 
time and refitting the 17-cluster scaling relations, shows 
that the cluster with the lowest ysph,5oo (Abell 383) af- 
fects the slope far more than any other. Without this 
cluster the rsoo scaling relation would have parameters 
A = 0.259;S5]3i, B = 0.56+°:!^, aM\Y= 20^^°%; in this 
case B is very consistent with the self-similar prediction. 
We have no strong reason to exclude this cluster, and the 
removal of outliers may bias the average scaling relation 
parameters, so no other outlier-excluded scaling relations 



are presented or discussed hereafter. The intrinsic scat- 
ter (TAf|F is found to be ~ 20% at all overdensity radii, 
which is significantly larger than the ^ 1 0% predicted in 
numerical simulations (e.g., Nagai|[2006 ) . This holds for 
the SS fits as well, confirming that the measured scatter 
is not strongly dependent on the slope parameter. 

The low scatter reported for simulated Y—M scaling 
relations does not capture additional complications intro- 
duced by the measurement techniques. One potentially 
important source of scatter in our data is the contribu- 
tion to l^ph from lower mass h aloes along the line of 



sight. This was examined by Holder et al. (20071 and 



found to be insignificant for our mass and redshift range 
The scatter between the weak-lcnsing-derived mass and 
true halo mass may also be important. In measurements 
of the type used here it i s expected to be ^ 20% (e.g., 
Becker & Kravtsov 2011), which is in good agreement 
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TABLE 3 

JW"wL — ^sph Scaling Relations 



Cluster 
Subset'^ 



B 



Cov(AB) 



^ (%) 



500 

1000 

2500 



All 
All 
All 



0.367"* 



7+0.096 
-0.099 
0^254+0.077 



0.118 



0.081 

+0.060 

0.066 



^+0.12 
-0.11 

a+0.11 
-0.11 

'^•"'^-0.13 



0.443 
0.48"* 



-0.012 
-0.009 
-0.008 



1-9 
-8 
^9 

20f? 



21"* 

l9^ 



500 


All 


n 9.1+0.036 
U.^4i_(, (,3g 


3/5 


500 


U 


n 907+0.043 


3/5 


500 


D 


n 14fi+0.062 


3/5 


1000 


All 


n 174+0.031 


3/5 


1000 


U 


212+0040 


3/5 


1000 


D 


104+0053 


3/5 


2500 


All 


ifi, +0.031 


3/5 


2500 


U 


n nfi+0043 

U.ilD_Q 343 


3/5 


2500 


D 


nfi9+0.062 


3/5 



20 
24 
17 
19 
24 
15 



TTtr 



13 

10 
+21 

9 
+9 

-7 

12 

9 

18 



19 
25 
22 



+9 

-7 

13 

9 

+24 

-11 



Note. — Scaling relations fitted according to Equation [6] For fits 
in the lower section of the table tlie power law slope was fixed to tire 
self-similar value (i?— 3/5). 

^ Cluster sample used in fitting. All clusters, (U) undisturbed, or (D) 
disturbed. 

The covariance of the parameters A and B. Because the scaling relation 
pivot point is not optimized, this quantity is non-zero and should be 
incorporated in any evaluation of the scaling relation uncertainty. 
^ Though the regression is performed in base-10 logarithmic space, these 
numbers are converted to base-e and are therefore percentage scatter. 

with the scatter we observe. 

5. DISCUSSION 

5.1. Comparison to Previous Results 

The M-wL— ^sph scaling relations presented above are 
the first to combine SZ measurements with weak gravi- 
tational lensing. Previous observational studies have as- 
sumed HSE to derive mass measurements, and conse- 
quently have different sets of systematics and biases. The 
assumption of HSE introduces an intrinsic correlation be- 
tween mass and Kjphj while the intrinsic correlation be- 
tween the lensing-based mass and l^ph in our study is 
expected to be lo wer (although not necessarily zero, as 
discussed in ^5.3 1. Moreover, systematic biases between 



HSE-based masses and lensing-b ased masses have been 

or limited at the 10% 
in cluster observations, 
masses and true masses 



measured ( Mahdavi et al. 



level dZhang et al.pOOH " 
and found between HSE- 



20081) 
UTDTi 



in simulated clusters (Rasia et a l. 2006; Lau et al.||2009 



Burns et al. 2010[ (Javaliere et al. 2011). I'hese Simula- 
tions suggest that nydrostatic masses are biased low by 
~ 10 — 20% at rsoo, with smaller biases at higher over- 
densities. In order to explore these effects, we compare 
our results with prev ious o bservational and theoretical 
studies in S 5.1.1 and ^5.1.2 respectively. 



5.1.1. Previous Observational Results 
Our results are consistent within ~ la with those of 



Bona mente et al.| (|2008|) at A = 2500 (left panel of Fig- 
[ip. |Bonamente et al.| analysed OVRO and BIMA 



ure 

observations of 19 clusters at 0.14 < z < 0.3 for which 
archival Chandra data were available. The clusters were 
modeled with an isothermal /?-model (after excising the 
central 100 kpc of the X-ray data) and masses were de- 
termined assuming HSE. The Y values were determined 
by normalizing the shape of the /3-model density (or pres- 



sure, since the cluster is assumed isothermal) profile with 
the SZ measurements and integrating. We calculate the 
Ygph values from the Markov chains used in that paper 
and perform the linear regression using the methods de- 
scribed above. Their results imply A = 0.131 ± 0.074, 
B = 0.65 ± 0.12, CTjvfiF = 9-5%7 which are consistent 
with self-similarity, and are statistically indistinguishable 
from our best fit parameters (Table Isl . Despite the lack 
of statistical signific ance, it is interes ting to note that the 
intrinsic scatter in Bonamente et al.^s HSE-based scaling 
relation is roughly half that found in our lensing-based 

relation. 

Our resul ts at A = 500 agree less well with | Andersson 



et al. ( 2011[ ) (right panel of Figure [T]) , who examined the 
-'Wsoo — ^sph scaling relati on for a sample of galaxy clusters 



discovered by the SPT ( [Vanderlinde et al. |2010| ). They 
estimated cluster masses from the pressure-like X-ray ob- 
servable Yx through an observational calibration of the 
-^^500 — Yx scaling relatio n with Mboq measured from 
X-ray data assuming HSE ( [Vikhhnin et al.||2009a| . The 
SZ profile, largely unresolved by the SP'l', was modeled 
using the observed X-ray density pro file and a universal 
temperature profile ( jVikhlinin et al. 2006). The result- 
ing scaling relation (A = 0.36 ± 0.03, B = Omtom, 
f^MlY = 9 ± 5%, in our nomenclature) is slightly dis- 
crepant with our FS model. However, comparison with 
our self-similar model is more direct, as Andersson et al. s 
slope matches the self-similar value. The scatter is again 
roughly half that observed in our data. Their normal- 
ization is ~ 30% higher in mass at fixed l^ph than our 
self-similar fit, a 2.3cr difference (after accounting for the 
contribution of the intrinsic scatter, compare the red and 
blue lines on the right panel of Fig IlJ . This difference is 
not attr ibut able to the hydrostatic mass bias discussed 
above (^5.1), because the assumption of HSE results in 
a bias in the direction opposite to the observed discrep- 
ancy. The origin of this difference remains unclear. 

5.1.2. Previous Simulations 

Accurate calibration of numerical simulations is a key 
step towards cosmological application of SZ surveys. 
Many authors have therefore developed simulations of in- 
creasing sophistication, five recent examples of which are 
over-plotted on our data in Figure [2] In all cases we show 
the predicted scaling relation between mass and Ygph, 
using fits provided by the authors for those simulation 
studies that did not publish a mass-ygph relation. The 
simulated relations are very close to self-similar {B = 0.6; 
Table |4]), and slightly under-predict l^ph at fixed mass 
relative to our best-fit self-similar model (red line in Fig- 
ure [2]). The fractional offset in l^ph of the simulated 
relations from the observed relation is ~ 10 — 35% at 
^^500 = 5 X 1O-^'*M0, the mean mass of the observed 
sample. We also rank the simulations using a simple x^ 
comparison of their predicted scaling relations with the 
observational data, assumingan intrinsic scatter of 20% 
in mass at fixed Ygph (Tahlem. We find that the 'Nag;iiJ 
( [2006; blue solid line in Figure [2]) simul ation is the clos- 
est ma tch to the observations, and the [Battaglia et al.| 
(20101 simulation is the poorest match, though not even 
~ 2(7 deviant from our self-similar relation. Note that 
this ranking does not depend on the adopted intrinsic 
scatter, although the absolute value of x^ does. 
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TABLE 4 
Simulated Mass— Kjph Scaling Relations 

A"" B^ Notes 



X 



2 c 



(1) 


0.27 


0.60 


(2) 


0.35 


0.55 


(3) 


0.32 


0.58 


(4) 


0.36 


0.54 


(5) 


0.37 


0.58 



M > IQi'^M, 



0. 



0,0.6 



M > 1.5 X IO^-^Mq, 0.1 < ^ < 0.2 
Preheating Model 
AGN Feedback Model 



15.8 
18.9 
17.8 
19.1 
26.0 



( l)|Nagai|l|2006); (2) |Sha 



et al.|l|2008l; ( 3)[Sehgal et al.| l |2010( ; 
(4) IStanek et al.| l |2010} ; (5) [Battaglia et ah] pOlO} . 

Scaling relations parameterized by A and B as in Equation Im 
^ X between model and the 18 data points (17 d.o.f) for an assumed 
intrinsic scatter of 20%. 

In a hybrid observational and simulated calibr ation 
of the A^wL— ^ph relation, Arnaud et al. (2010) con- 
structed an "observed" scaling relation between Ygph and 
-Afsoo without SZ measurements. They use the univer- 



sal pressure profile described in Section [3T2| which they 
calibrate against X-ray pressure profiles and simulated 
clusters, inside and outside of rsoo, respectively. Their 
best-fit relation (^=0.34, B=0.56) also app ears in Fig- 
ure [H It m atches well with the s imulations of Sha w et al.l 
(2008) and|Sehgal et al. (2010), similarly underpredict- 
ing y at fixed M. 'I'he mass scale for this scaling rela- 
tion comes from their M—Yx scaling relation, which is 
calibrated against hydrostatic masses, so the line might 
instead have been expected to over-predict Y due to the 
hydrostatic bias. At the measurement precision and level 
of scatter in this work we do not detect such a bias, even 
though it is expected to be more important at rsoo than 
at r'2500 where we are ab le to compare with the results of 
Bonamente et al. ( 2008 ) . 



5.2. Observed Scatter and Morphological Segregation 

The distribution of fractional deviations in mass 
(AMa/Ma) of clusters from the mean Ma — l^ph re- 
lations is asymmetric, particularly at A = 500 (Fig- 
ure pi) , where it peaks to the low-mass side of the re- 
lations (negative deviations) and has an extended tail to 
the high-mass side (positive deviations). Interestingly, 
the clusters classified as "disturbed" on the basis of the 



centroid shift measurements discussed in [3.3 all have 
^Afsoo/Msoo < 0, and the tail of the distribution at 
^Afsoo/Msoo > 0.3 comprises solely undisturbed clusters 
(compare the blue and red in Figure |3]) . This suggests 
that the large intrinsic scatter in our mass-Ksph relations 
may be related to cluster morphology. 

We detect a significant difference in the normalization 
of the self-similar mass-Vsph relations between disturbed 
and undisturbed clusters at all three overdensity radii. 
At fixed Ygph, the mass of undisturbed clusters exceeds 
that of disturbed clusters by 13 ± 6%, 28 ± 5%, and 41 ± 
6% at A = 2500, 1000, and 500, respectively (Table [s) 
Figure |4| . The probability of observing such an offset 
is low even in scattered data like these; drawing random 
samples (with replacement) of 6 and 12 clusters from 
our data and repeating our fits, we find that only 1% 
of random sub-samples are more significantly offset than 
the real data at r5oo. 

If confirmed as a real physical effect, the morphologi- 
cal segregation of clusters in the mass-Fsph plane would 
have major implications for SZ surveys, which expect to 
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Fig. 2. — Data and scaling relations at A = 500 compared to 
scaling relations measured from X-ray data and numerical simu- 
lations of galaxy clusters. The scaling relation fits, plotted as in 
Figure IT] are truncated at the boundaries of the mass range of 
the cluster sample including Icr uncertainty. Data points from this 
work are shown as open circles without error bars. Five simulated 
scaling relations between Vsph and true M500 arc shown, as well as 
the scaling relation predicted by [Arnaud et al. (20101 from X-ray 
observations and numerical simulations, ihe INagai (2006) scaling 
relation was derived fr om simulate d clusters with mass greater than 
10^* Mq, while the [Sehgal et al.|(2010J relation was determined 



from simulated M5oo>2xlCF 
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Fig. 3. — Distribution of fractional deviation in mass at A = 500 
of the 18 clusters in our sample from the best-fit self-similar scaling 
relation fit shown as the solid red line in th e right panel of Fig. n] 
Disturbed clusters (blue; defined in j ]3.3| l lie exclusively on the 
low-mass side of the mean relation; uridislurbed clusters (red) lie 
on both sides of the mean relation, with a tail to large positive 
deviations. 
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Fig. 4. — The mass-y scaling relation at each overdensity, with cluster morphology indicated by the data symbols. Disturbed clusters are 
marked with open squares and undisturbed clusters with filled circles. The self-similar fit for the full sample is shown with a solid (black) 
line, and the self-similar fits to the undisturbed (red dashed) and disturbed (blue dot-dashed) subsamples are also shown. 



produce nearly mass-limited cluster samples minimally 
biased with respect to dynamical state. However, mor- 
phological segregation has not been seen in mass-ygph 
relations predicted from numerical simulations. Indeed, 
it has generally been found that the SZ signal is a ro- 
bust mass proxy (merging and non-merging clusters are 
co-located in the ma ss-ygph plane) even during periods o f 
mass accretion (e.g., Motl et aL||2005 Poole et aL||2007 ). 
Simulations of merging clusters do hnd that cluster 
mergers produce transitory boosts in the SZ signal, al- 
though when integrated over t he entire cluster these 
boosts are small and short -lived (IPoole et al. | |2007 Wik 
et al.||2008). For example, Poole et al. showed that tor 



most merger mass ratios and relative velocities, the ob- 
served Y is usually below the final "post-merger" value 
during the merger. This is a natural consequence of the 
finite time required to heat the ICM to the higher equi- 
librated temperature of the post-merger halo. The ob- 
servational signature of this effect is that at fixed post- 
merger halo mass, merging (disturbed) clusters should 
have Y suppressed relative to non-merging (undisturbed) 
clusters. These theoretical predictions are in the oppo- 
site sense to the morphological segregation detected in 
our observed mass-l^ph relations. 



5.3. Correlation of Mass and Y Measurements 

The ideal scaling relation measurement compares two 
observables that have been derived independently. X- 
ray scaling relations have typically been constructed from 
HSE-based masses derived from the same X-ray data that 
appear in the observable axis. HSE-based scaling rela- 
tions therefore inevitably suffer varying degrees of cor- 
relation that depend on the details of the measurement 
methods and the observable a gainst which mass is plot- 



2010a|. This correlation acts to 
To 



erred ior HSE-based scaling rela- 



ted (e.g., Mantz et al. 
suppress the scatter int' 
tions. 

The MwL and Vsph measurements presented in ^ 
are also not completely independent, as the outer in- 
tegration boundaries for the SZ profile are set by the 
lensing-derived overdensity radii. This introduces a cor- 
relation between observables. The mass error (SM = 
MwL — Mtrue) for an individual cluster corresponds to 



an error in the overdensity radius (Sr), 

SM = 47rr^p(r)(5r , (7) 

and this error is transferred to the calculation of ygph- 
Equation [7] can be converted to an equation involving 
fractional errors in M and r by dividing by M{r) = 
47rrV/3, 

d{log{M)) = 3^6{log{r)) . (8) 

The perturbation in the calculated Ygph introduced by Sr 
follows the same form as equation [8| with the pressure 
substituted for p and l^ph for M. 

From the equations for the logarithmic radial slopes 
of M and Ksph we can infer the response of these two 
quantities to the errors in ta- The direction of motion 
in the Mwl— ^sph plane for an error Sr, expressed as the 
slope Slog{M) / SlogiY) , is just the ratio of the logarith- 
mic slopes in density {Slog(M) / Slog{r)) and pressure at 
the radius of interest. This ratio varies from cluster to 
cluster and across overdensities, in the present sample, it 
is (on average) 2.3 at ^2500, but 1.7 at r^QQ. The latter is 
very close to the slope of the self-similar scaling relation, 
1/B = 5/3. On average then, the scatter between Mwl 
and Mtrue leads to errors in the derived l^ph that move 
the clusters along the scaling relation at r^QQ. 

It was noted in §4.3| that the scatter between the mass 
derived from the WL shear profile fitting techni que and 
the true cliister m ass is expected to be ^-^20% (jBecker 
fc Kravtsov| 2011 1, which matches the values in TaBIeW 
I'hat the observed and predicted scatter agree so well 
may be coincidence, simulations suggest that the astro- 
physics of the ICM can contribute another 10-15% scat- 
ter between Y and true mass, and systematic effects like 
those discussed in subsequent sections can be expected 
to further increase the observed scatter. The above anal- 
ysis suggests that aM\Y could be larger if the Ygph and 
Mwl measurements were decoupled. The clusters in this 
sample all have Chandra data, which provides us with an 
independent measurement of rsoo for the calculation of 
yjDh. Using the X-ray values of rsoo from by Sanderson 



et al. ( 2009 ) as integration radii for both the iNl^'W-lialo 
fits to the weak lensing data and the SZ profiles, we again 
find the 20% scatter in the scaling relation. When these 
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Fig. 5. — X-ray centroid shift versus slope of the gas density pro- 
file at 0.04rK nn (|o|)- The density slopes are taken from Sanderson 
|et al.| l [2009[ l, where a < —0.85 provided a reliable dividing line 
between cool-core and non-cool core clusters. 

radii are used for the SZ data alone, the scatter increases 
to 26%. The couphng between observables through the 
integration radius therefore does appear to lessen the ob- 
served scatter, and the scatter of fully independent mea- 
surements may be more representative of the true intrin- 
sic scatter. 

5.4. Choice of Pressure Profile 

Reliable measurements of Ygph depend on the appropri- 
ate choice of pressure profile shape in the ICM models. 
Systematic errors in the choice of pressure profile may 
propagate to systematic errors in Ygph measurements, 
and thus contribute to the observed scatte r and segrega- 
tion. Our l^ph measurements employ the Arnaud et al. 
(2010) average pressure profile. However, these authors 
noted that CC clusters have more concentrated pressure 
profiles than disturbed clusterfP^ 

As a test of the importance oTthe profile choice on our 
SZ measurements, we recalculate our Ysp^ values u sing 
the morphology-s pecific pressure profiles presented in |Ar-| 
naudetal. (2010). For cluster s we classify as "disturbed" 
and "undisturbed" we use the |Arnaud et al.| "morpholog- 
ically disturbed" and "cool-core" profiles, respectively. 
The change in profile leads to small systematic changes, 
the disturbed clusters decrease in Ygph by 1% and 3% 
at ''2500 and ^500, respectively, and undisturbed clusters 
increase by 1% and 5% at these radii. These changes 
reduce the morphological segregation in the M500 — i^sph 
relation from 41% to 31%. 

We expect this decrease in morphological segregation 
to be an upper limit on the effect of pressure profile choice 
on our results. Although we use the Arnaud et al. (20101 
pressure profile for all of our undisturbed clusters, not 

^^ This cool-core/disturbed classification identifies "disturbed" 
clusters on the basi s of (w), and CC clusters on the basis of their 
central gas density dPratt et al.|2009||. 



all of them host a cool core (Figur e Isl). In the o rigi- 
nal measurement of the CC profile, Arnaud et al. ex- 
cluded undisturbed clusters that lacked cool cores from 
their subsample, and these exclud ed c lusters had profiles 
closest to the one we adopted in §3.2| We therefore con- 
clude that although up to one foiirth of the morphologi- 
cal segregation can be attributed to systematic errors in 
Ysph due to incorrect pressure profile choice in the ICM 
models, the majority of the segregation remains to be 
explained. 

5.5. Choice of Weak-lensing Mass Model 

Reliable measurements of Mwl also depend on appro- 
priate model choice. The weak-l ensing masses used in 
this article are based on the initial [OklO| analysis of data 
from our Subaru observing program, in which the mass 
distribution of each cluster was descri bed b y a model con- 
taining a single spherical NFW halo ([ 3.1 1. We therefore 
explore whether this choice may contribute to the mor- 
phological segregation in the Mwl— ^spii plane. 

Additional halos that might be justified by the data in- 
clude sub-halos within each cluster, and large-scale struc- 
ture projected along the same line-of-sight. We attempt 
the lowest order correction to our lensing analysis by 
adding a second dark matter halo to each mass model 
at the position of the most significant sub-peak i n each 
of the density maps presented in Appendix 1 of |OklO 
These 2-halo models were fitted to the 2-dimeiisiona l 
weak-shear field using lenstooiJ^ (Julio et al. 2007). 
Comparing these models to the original 1-halo mod- 
els, we found that the Bayesian evidence ratio signifi- 
cantly favored a two-halo model for all but 3 clusters. 
We re-calculated AfwL at A = 500 for each cluster us- 
ing the more probable of the two models, excluding the 
sub- halo if it is known to lie outside the cluster (e.g., 
RXJ0153.2-h0102 adjacent to A 267; see Figure 21 in 
OklO ), and assuming that the projected separation of 
the two halos is the true three-dimensional separation of 
the halos in all other cases. In the resulting M^QQ-Ygph 
scaling relation we find that neither the scatter nor seg- 
regation are reduced within the uncertainties. 

The one-dimensional radial profiles used to model the 
clus ter weak-lensing si gnal may also be a source of scat- 



ter. [Oguriet al.| ( |2010[ ) fitted elliptical weak-lensing mass 
models to UklO s data, thus relaxing the assumption of 
azimuthal symmetry in the single NFW halo of the lat- 
ter's models. We find that the scatter and segregation 
in the MwL~ ysph relation a re all unchanged within the 
uncert ainties if |Oguri et al.[ s masses are substituted for 
OklO s masses. Both of these tests suggest that the de- 
scription of the projected mass distribution is not the 
dominant source of scatter. 

The line-of-sight structure of the cluster dark mat- 
ter is predicted to introduce systematic errors to lensing 
mass es derived assuming a spherically symmetric p rofile 

iMeneghetti et al ||2010| . An 



(e.g., Corless fc King||2007 



observational indicator of t 



le line-of-sight elongation of 



the mass distribution would therefore be a powerful tool. 
First, we consider the X-ray centroid shift measurements 
upon which the disturbed/undisturbed classification was 
made in j ]3.3[ We plot centroid shift versus fractional 
deviation in mass from the best-fit SS scaling relation in 



■^^ http : //www . oamp . f r/cosmology/lenstool/ 
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Fig. 6. — Fractional deviation in mass from the SS Afsoo-K^ph 
relation versus centroid shift {w) (left) and ellipticity of the BCG 
in each cluster (right). Filled circles (red) and open squares (blue) 
represent undisturbed and disturbed clusters, respectively. The 
classification is done according to the centroid shift value and the 
dividing line is marked by the vertical dashed line in the left panel. 

Figure [6J finding no reliable trend between the mass er- 
ror and (it;). Though the large-(u;) clusters (disturbed) 
lie exclus ively on the low-mass side of the relation, as 



noted in ^5.2 the small- (w) (undisturbed) clusters span 
a wide range in deviation (AAfsoo/Msoo '^ — 0.75) 
at all centroid shift values. The largest deviations sug- 
gest that such clusters may be prolate s pheroids viewed 
along an axis close to their major-axis dCorless & King 



20071, while the smallest deviations suggest that the 
dark matter in such clusters is relatively undisturbed and 
spherical. This implies that the interpretation of small 
centroid shifts is degenerate between a relatively undis- 
turbed ICM and the disturbed ICM of a prolate cluster 
whose major axis is closely aligned with the line-of-sight. 
A "cleaner" indicator of the orientation of cluster dark 
matter halos is therefore required. 

BCGs are prolate stellar systems with their major axis 
closely aligned with the major axis of the cluster dark 
matter halo and ICM (e.g., Hashimoto et al. 2008 Fasano 
et al. 2010). The major axis of a prolate BCCJ that is 
measured to be strongly elliptical in projection on the sky 
likely lies close to the plane of the sky, and thus so too 
does the major axis of the cluster dark matter halo. An 
almost circular BCG, on the other hand, indicates that 
the major axis is close to the line-of-sight through the 
cluster. BCG ellipticity measurements are available from 
surface photometry of Hubble Space Telescope (HS T) ob- 
servations o f 15 of the 18 clusters in our sample (Smith 



et al. 



2010). We analyzed two of the remaining clusters 
0142.0-1-2131 and A 697) in the same manner (the 
ZwCl 1459.4+4240, has not been observed 



last cluster, 

with HST). We find that the clusters with the roundest 
BCGs suffer the strongest positive deviations from the 
-^^WL~^ph relation, and that clusters with the most el- 
liptical BCGs suffer negative deviations (right panel of 
Figure [6|. Most significantly, clusters with ebcg < 0.15 
have exclusively large deviations (AM500/M500 ^ 0.5). 
The mean deviation of these "round" BCGs (e < 0.15) 



is (AM500/M500) = 0.76 ± 0.17 and of "elliptical" BCGs 
(e > 0.15) is (AM500/M500) = -0.10 ± 0.17, where the 
error bars are the la standard deviations on the means. 
The correlation between AM500/M500 and eecG is 
strikingly similar to the trend between Mao wl/Msd true 
and viewing angle in the numerically sim ulated clus- 
ters presented in [Mcneghet ti et al.| (2010). Indeed, vi- 
f pe 



sual inspection ofplencghetti et al.[s I'lgure 17 suggests 



that we are viewing clusters with "round" BCGs within 
^ 20 — 30° of the major axis of the cluster inertia ellip- 
soid, and with "elliptical" BCGs at larger viewing angles. 
An observed BCG ellipticity of eecG = 0.15 (correspond- 
ing to an observed axis ratio oi q = b/a = 0.85) can be 
converted to a viewing angle with respect to the major 
axis by assuming an intrinsic BCG axis ratio ( /3). We 
adopt J3 = 0.67, following the measurements of [Fasano] 



et al. (2010). The viewing angle is then estimated to be: 



Tp = arccos 




w/qy 



/32 



34° 



(9) 



which is consistent with the visual compa r ison o f our 
Figure [6] and Figure 17 of |Meneghetti et al.| ( [20T0 |. 

Our oata provide observational support for previous 
numerical studies that have highlighted projection effects 
as an important source of systematic uncertainties in 
weak-lensing cluster mass measurements. These uncer- 
tainties arise from both the underlying shape of the main 
cluster halo and sub-halos that reside within the cluster 
due to a merger and/or in the surrounding large-scale 
structure. Future progress on the use of weak-lensing to 
calibrate mass-observable scaling relations will therefore 
benefit from intrinsically asymmetric lens models that 
include halo triaxiality and/or multiple halos. Direct 
measurements of the ratio of the line-of-sight depth to 
the plane-of-sky size of the cluster through joint X-ray 
and SZ measurements of the ICM may aid in construct- 
ing a more three-dimensional halo model. These efforts 
should lead to a reduction in the scatter and segregation 
in the A/wl — i^ph relation. 

6. CONCLUSIONS 

Surveys of galaxy clusters using the SZ effect require 
a well-calibrated mass-observable scaling relation in or- 
der to constrain cosmological parameters. In this paper 
we report the first comparison of the SZ observable Ygph 
with the weak lensing-determined halo mass. We find a 
strong correlation, and our best-fit scaling relation is con- 
sistent with expectations from self-similarity arguments 
(Table [3]) . Comparing to previous scaling relations mea- 
sured against hydrostatic mas ses, our normalization at 
^"2500 is consistent with that of Bonamente et al.| (|2008|) 
bu t at rson it mildly differ s from SPl' results presented 
by Andersson et al. (2011). The scatter in our relation 
is larger than predicted in simulations of the SZ signal, 
though it is not unexpected given the difficulties inher- 
ent in deriving spherical masses from the natively two- 
dimensional weak lensing measurements. 

We considered the morphology dependence of the 
A^WL— ^sph relation by dividing our sample into undis- 
turbed and disturbed sub-samples based on the well- 
established centroid shift parameter. We found that 
clusters segregate in the Mwl— i^ph plane: the mass of 
undisturbed clusters exceeds that of disturbed clusters 
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at fixed i^sph- This effect becomes more pronounced at 
larger radii (see Table |3|, at A = 500 undisturbed clus- 
ters are 40% more massive than disturbed clusters. Such 
a segregation is not predicted by numerical simulations 
for comparisons of the SZ signal and true cluster mass. 
We discuss a wide-range of possible interpretations of the 
segregation, and conclude that the relative simplicity of 
the models fitted to the data is the most likely cause. No 
more than one fourth of the segregation can be attributed 
to morphological biases introduced by the adoption of a 
single (morphologically blind) pressure profile in the SZ 
data analysis. From multiple lines of evidence we con- 
clude that most of the segregation is caused by modeling 
intrinsically prolate cluster dark matter halos as spheri- 
cal objects when determining weak lensing masses from 
our data. Future uses of weak lensing as a calibration 
tool will be aided by attempts to correct for halo ori- 
entation, and we will explore this in a larger sample of 
objects in a future paper. 
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